o 
o 

(N 
> 

o 

(N 



o . 

"—I . 
-(— > ' 

Oh' 

q: 

o ■ 

•I— I ■ 

Oh, 



> 

in 

(N 



O 



X 



Unified Homogenization Theory for Magnetoinductive and Electromagnetic Waves in 

Split Ring Metamaterials 
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A unified homogenization procedure for split ring metamaterials taking into account time and 
spatial dispersion is introduced. The procedure is based on two coupled systems of equations. The 
first one comes from an approximation of the metamaterial as a cubic arrangement of coupled LC 
circuits, giving the relation between currents and local magnetic field. The second equation comes 
from macroscopic Maxwell equations, and gives the relation between the macroscopic magnetic 
field and the average magnetization of the metamaterial. It is shown that electromagnetic and 
magnetoinductive waves propagating in the metamaterial are obtained from this analysis. Therefore, 
the proposed time and spatially dispersive permeability accounts for the characterization of the 
complete spectrum of waves of the metamaterial. Finally, it is shown that the proposed theory is 
in good quantitative and qualitative agreement with full wave simulations. 



PACS numbers: 

Diamagnetic properties of systems of conducting rings 
are known since long by physicists. In 1852 Wilhem We- 
ber [ij tried to explain natural diamagnetism (discovered 
by Faraday some years before) as a consequence of the 
excitation of induced currents in some hypothetical con- 
ducting loops that supposedly existed in diamagnetic ma- 
terials. In order to enhance the magnetic properties of ar- 
tificial media (or metamaterials in modern terminology) 
made from metallic conducting rings, S. A. ShelkunofF 
proposed in 1952 to introduce a capacitor [2|, so as the 
rings become resonant. More recently J. B. Pendry et al. 
Q proposed to replace the capacitively loaded rings by 
planar split ring resonators (SRRs) which substitute the 
lumped capacitor by a distributed capacitance between 
the rings. Because Pendry's SRRs can be easily manu- 
factured by using standard printed circuit technologies, 
this design opened the way to manufacturing true mag- 
netic metamaterials made of many individual elements 
(SRRs) at many laboratories around the world. As a 
consequence of its resonant behavior, capacitively loaded 
rings and/or SRRs can produce metamaterials with neg- 
ative magneticpermeability above resonance. It is also 
well known [J, |5|, |6|] that when a system of these elements 
is properly combined with another system of elements 
(metallic wires or plates, for instance) producing a nega- 
tive electric permittivity [7|, a metamaterial with simul- 
taneously negative permittivity and permeability (or left- 
handed metamaterial [8|) arises in the frequency band 
were both subsystems present negative parameters. Re- 
markably, the electric and magnetic properties of such 
combinations are, quite approximately, the superposition 
of the electric and magnetic properties of each subsys- 
tem. This superposition hypothesis is not apparent at 
all, since the elements of both subsystems must be placed 
closely in space, and therefore electromagnetic coupling 
may be present. In fact, this has been a controversial 



issue (see for instance [3, [ifl]). For the specific SRRs 
and wires configuration proposed in [J] this hypothesis 
has been recently demonstrated by one of the authors 



ill ] and, in general, it can be admitted that the afore- 



mentioned superposition hypothesis is valid provided the 
elements of both subsystems are placed in such a way 
that their quasi-static fields do not interact or interact 
weakly [lO|. Almost simultaneously, other analyses and 
experiments 1^ ll3| did show that SRR based metama- 
terials also support, in some frequency bands, slow waves 
based on short range interactions between the SRRs; the 
so called magnetoinductive (MI) waves, which can not be 
deduced from the usually assumed local magnetic perme- 
ability of the metamaterial. Interestingly, many of the 
physical effects expected in negative permeability and 
left-handed metamaterials, such as frequency band gaps 
and frequency bands of backward- wave propagation, also 
come out from the analysis when the coupling between 
electromagnetic and MI waves in SRR systems is consid- 
ered 1J|, thus providing an alternative exp lanation for 
such effects. Although the analysis in [lj| has a great 
heuristic value, it can not be considered as fully satisfac- 
tory because it only considers one-dimensional systems in 
the nearest neighbors approximation. On the other hand, 
the presence of waves which can not be deduced from a 
local time dispersive magnetic permeability in split ring 
metamaterials can be expected from the fact that its pe- 
riodicity is usually not smaller than one tenth of a wave- 
length. As it is well known [15| , when the periodicity of a 
given medium approaches the wavelength of the electro- 
magnetic radiation, it becomes not only time dispersive 
but also spatially dispersive. Therefore, it can be ex- 
pected that both, electromagnetic and MI waves, would 
come out from the analysis if spatial dispersion in split 
ring metamaterials were taken into account. In fact, the 
main purpose of this paper is to show that both kinds 




(a) 



(b) 



FIG. 1: (a): Unit cell of a material formed by a cubic array 
of current loops. Each unit cell has three current loops cen- 
tered at the faces of the cube, (b): Similar to (a) but for a 
realistic metamaterial formed by edge coupled SRRs with two 
splits. The SRRs are formed by circular wires with radius r. 
The distance between the inner and outer rings is d and the 
average radius of the particle is R. 



equations lead to 



e) Ho + (fc,2„ - kk.) Mo = 0, 



(3) 



with fc,„ = koy/e^. = ujy/erSQfio, where Sr is the macro- 
scopic relative dielectric constant of the metamaterial. 
By combining ([T|) and ^ the following equation for lo is 
obtained 



Z(k, w) • lo = -juA^iQ Ho 



Mo 



(4) 



where Z(k, uj) is an impedance matrix which incorporates 
all the magneto-inductive effects between the neighboring 
rings. Explicit expressions for the diagonal and the off- 
diagonal terms of Z are 



of waves can be obtained from a unified analysis when 
a proper homogenization procedure, taking into account 
spatial dispersion, is developed. 

In order to simplify the analysis, we will consider an 
ideal metamaterial made of a cubic arrangement of LC 
circuits supporting current loops as sketched in Fig. la 
for a unit cell. We define the current vector I" on each 
unit cell as I" — {I^,Iy,I^), where n = {nx,ny,nz) 
specifies the location of each unit cell in the lattice and 
/" denotes the current along the loop located in the face 
normal to the i-direction of unit cell of index n. The time 
dependence is assumed of the form I" ex exp{JLut). Each 
component of the current vector is governed by equation 



jujL 



1 



jujC 



/" 



-JUJ'^?, 



(1) 



where L, C are self-inductance and self-capacitance of 
the circuit (losses are neglected by simplicity, although 
they can be easily introduced in the analysit through a 
ring resistance) . In Eq. ([1]) (f>f is the total magnetic fiux 
through the considered loop which, using Lorentz local 
field approximation, can be calculated as 



$r = ^Mo H, 






r<R 

E 

(m#n)A(j#i) 



*J 3 



(2) 



where A is the area of the loop, H, M are the macro- 
scopic magnetic field and magnetization on the ring, and 
M"™ are the mutual inductances between the loops with 
indexes n and m. The summation is extended only to a 
sphere around the unit cell n which contains several unit 
cells but which radius R is smaller than the wavelength. 
In the following we will assume a spatial field depen- 
dence of the kind {H, M} — {Hg, Mq} exp(— jk • r) and 
I" = Iocxp(— jak • n), where a is the lattice periodicity 
and lo. Ho, Mq are constant vectors. With the assumed 
time and space dependence, the macroscopic Maxwell 



Zu = jujL < 1 



Wq 



n/0 



-jakn 



(5) 



z,. 



Zji 



r<R 71 rOn 



* ^ J (6) 



where ojq = l/vLC is the frequency of resonance of the 
rings, with similar expressions for the remaining compo- 
nents of Z- Taking into account that Mo = AIq/q^, it is 
possible to combine ^ with (j4]) which gives 



ziKuj) 



jojUoA^ 2k'i + P- 3kk 1 
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• In = 0. 



(7) 



The dispersion equation for plane waves in the metama- 
terial can be obtained by equating the determinant of 
([T]) to zero. In the most general case this equation can 
be only solved numerically. However, in some cases, it 
is possible to give analytical solutions. In particular, for 
propagation along one of the coordinate axis (for k = fcic, 
for instance) the summation in ([6|) vanishes because all 
unit cells in planes perpendicular to the a;— axis are in 
phase and mutual inductances cancel out couple by cou- 
ple (for instance, the mutual inductances between the 
ring marked Ix , and the rings placed on the top and lower 
faces of the cube of Fig.l.(a) cancel each other, and so 

on). Therefore, the matrix Z(k, cj) becomes diagonal, 
and ([7]) can be easily solved. This gives two branches: a 
longitudinal wave with lo = Iqx^ given by 



2 jujfioA _ 

^xx ^ — ^ — U 



3 a' 
and a transverse wave with lo 

JLuiipA^ 2kl + kl 

Otl '^'tri "^T 



3 - -^ (8) 

= hyY + lozZ given by 



0. 



(9) 



If only interactions with the closest rings are considered 
for the computation of the summations in ^ and ([6]) the 
dispersion equation for the longitudinal wave becomes 



= 1 



2—— cos(aA:a;) 



Mr 



(10) 



2ao 

L 3a3' 

where ao = IJ^q^ ji-i and Ma and Mc are the mutual in- 
ductances between closest rings of the same orientation, 
placed in the coplanar and the axial directions, respec- 
tively (as it was already mentioned, the contributions of 
the mutual inductances between rings of different orienta- 
tions cancel each other for this particular propagation). 
It can be easily recognized in (|10p the dispersion rela- 
tion for longitudinal magnetoinductive waves [12l |. with 
some small corrections, which take into account the ef- 
fects of the rings other than the nearest neighbors in the 
axial direction. In the same approximation the dispersion 
equation for transverse waves can be written as 
1,2 Ha 

1.2 



1 



2(A/a+Afe) 



_Q0 



For high values of fc. 
2M. 



; ^ (11) 

(fcj. 3> fc„i) this equation reduces to 



-^ = 1 + 

9 



cos{akx) 



2{Ma + M,) 



V-^-^ (12) 



which corresponds to the dispersion relation for trans- 
verse magnetoinductive waves [12l |. On the other hand, 
in the long wavelength limit [akx <C 1) (fTTj) reduces to 

jU2 



x(^) 



1 = 
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2Af„ 
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a a 



(13) 



This equation gives the value for the magnetic suscep- 
tibility that is obtained when the Lorentz homogcniza- 
tion procedure is applied to the metamatcrial with the 
ring magnetic polarizabilities a = aQiuj'^/uj^ — 1)^^ al- 
ready proposed in [l6| , except for a small correction term, 
2Ma/L + 4Mc/L, accounting for the effect of the closest 
rings. Actually, if such correction term is calculated by 
assuming a magnetic dipole approximation for the rings, 
it can be easily shown that it vanishes, thus giving ex- 
actly the Clausius-Mossotti formula for the susceptibility. 
Therefore, we can conclude that in the long wavelength 
limit, the transverse waves © correspond to the electro- 
magnetic waves that are obtained from the local time- 
dispersive permittivity /i — ^o{l + x{<^)}- Conversely, 
in the short wavelength limit {k^ 3> km), they converge 
to the transverse magnetoinductive waves (fT2|) . Further- 
more, from Mq — Alo/a^ and Q equation 



Mo 



f Ja^ =., N 1 



Hn 



X(k,c^).Ho (14) 



can be obtained. Now in ^ Mo can be replaced by H 
leading to 



{(■ 



e 



-kk + ki 



X(k,c^)} 



Ho = 0, (15) 
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FIG. 2: (Color online) Dispersion diagram along F-X, F-M, 
F-R directions obtained from the analytical model (lines 1, 2, 
3) and from a full- wave simulation (dashed line). The coor- 
dinates of selected points are F = (0,0,0), X — (7r/a,0, 0), 
Af — {ir/a, 7r/a, 0), R — (vr/a, vr/a, n/a). 




FIG. 3; Dispersion diagram along the closed path R-F-X-M-R 
obtained from full-wave simulation. Horizontal axis shows a 
projection of k on the corresponding line on the boundary of 
Brillouin zone. 



which gives the same dispersion equation as ([7]). There- 
fore we can conclude that the non-local (i.e. time-and 
spatially-dispersive) magnetic permeability 



M(k,cj)=^o 1 + 



-,3 _ 



fioLuA-' 3 



(16) 



provides a complete characterization of the metamatcrial, 
accounting for all kind of waves propagating through it. 
In the long wavelength limit (a|k| ^ 1) all the exponen- 
tial terms in Z(k, lu) can be equated to unity and, from 
the aforementioned cancellation of the inductances be- 
tween rongs of different orientation, all the off-diagonal 
terms of Z(k, w) vanish. Thus, Z(k, uj) becomes an scalar 
and the magnetic susceptibility in (fH)) reduces to the 
scalar x(a;) in ([13]). 

As a numerical example, we have studied the propaga- 
tion of electromagnetic waves in a metamatcrial formed 
by the simple cubic lattice of split ring resonators whose 
unit cell is depicted in Fig. lb. The capacitance C and 
self inductance L were calculated following the ideas of 



Ref. lr| but for the case of an SRR made of wires instead 
of planar strips [22| • The mutual inductances M"™ were 
calculated numerically using Neumann's formula includ- 
ing time retardation. The macroscopic permittivity was 
approximated by substituting the SRRs by planar con- 
ducting disks of the same external radius. Permittivity 
of such medium was calculated using static Lorentz ho- 
mogenization theory [l8( and its value is Er — 2.5. Using 
Eq. ([7|) and the first neighbors approximation, the dis- 
persion characteristic of the electromagnetic modes sup- 
ported by the metamaterial along different directions of 
the first Brillouin zone were calculated. The result for 
the geometry associated with R = 0.44a, r = 0.005a and 
d = 0.03a (see Fig. lb) is depicted in Fig. 2. It can be 
seen that the band structure is formed by three branches 
and contains a bandgap for all the depicted directions of 
k. For or k along F-X the first and third branches cor- 
responds to the transversal mode described by Eq. ([9]) 
and the second branch corresponds to the longitudinal 
mode described by Eq. ([5]). Figure 2 also shows the 
high isotropy of the transversal mode even for moderate 
values of k, a fact that is expected from the tetrahedral 
symmetry of the system |19l| . 

To assess the accuracy of the proposed analytical 
model, we have also numerically computed the exact 
band structure of the aforementioned periodic material 
using the hybrid-plane- wave-integral-equation formalism 
introduced in I2u\ . The result of the numerical simulation 
is presented in Fig. 2 for the specific direction F-R, and in 
Fig. 3 for the closed path R-F-X-M-R. Good qualitative 
agreement between theory and simulation can be seen 
from Fig. 2. We think that the quantitative disagreement 
for high values of k in the second branch can be attributed 
to the specific local field approximation considered in Eq. 
2, which is strictly valid only for small values of k. In the 
case of the third branch, the disagreement is due to the 
proximity of the second resonance of the SRRs, which 
is not taken into account in the model. This effect is 
more visible in Fig. 3, where higher frequency branches 
are included. This figure also shows a complete elec- 
tromagnetic band gap in the range 1.18 < kga < 1.50, 
in agreement with the hypothesis of a negative perme- 
ability in such frequency band. It may be worth noting 
that the effect of the substitution /c^ -^ ~km (or equiv- 
alently e,. -^ —£r) into (jlip is the onset of a backward 
wave passband in the frequency range of the stop-band 
of Fig. 3, as well as the conversion of the passbands of 
Fig. 2 into stop-bands. Therefore, it can be guessed that 
the proposed model will be also useful for the analysis 
of isotropic left-handed media made of SRRs and wires 
or any other elements providing a macroscopic negative 
permittivity (provided the conditions for the validity of 
the superposition hypothesis previously discussed in [10| 
are fulfilled). Work in this direction is in progress. 
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In this case the total capacitance of the SRR can be esti- 
mated as C = ^nRCpui = Tv^eoR/ (4cosh"^ f ^ ~ -'^))' 
where Cpui is the per unit length capacitance of two par- 
allel wires given in ^l|(p. 88). The self-inductance of 
the SRR was approximated by the self-inductance of a 
torus with large radius R and small radius r [2j|(p. 233): 
L = fioR [In (8f ) - 2] . 



